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Abstract 

The exact solution is obtained for the eigenvalues and eigenvectors of the 
row-to-row transfer matrix of a two-dimensional vertex model with unlimited 
number of states per bond. This model is a classical counterpart of a quantum 
spin chain with an unlimited value of spin. This quantum chain is studied us- 
ing general predictions of conformal field theory. The long-distance behaviour 
of some ground-state correlation functions is derived from a finite-size analysis 
of the gapless excitations. 
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Most of exactly solved models in two dimensions are known to be related to the six- vertex 
model solved by Lieb |IJ. This model exhibits an infinite number of conserved integrals 
commuting with each other. The generator of these integrals of motion is the row-to-row 
transfer matrix. The conserved quantities are quantum Hamiltonians describing different 
one-dimensional systems, in particular one of these describes the XXZ chain [2]. There 
exist also exact solutions for a variety of generalizations of this model. In particular lattice 
models where the link variables may assume q > 2 different states were solved using the 
Bethe-anstaz method [3-12]. 

In this letter we present an exactly integrable vertex model in which the number of 
different states of a given link is not limited, and may be considered as a generalization of 
the six-vertex model for the case where on each link of the lattice we may have an unlimited 
number of lines. This model is the classical counterpart of a quantum spin chain with 
unlimited value of the spin. The vertex model may be formulated as follows. We consider 
a square lattice with M rows and iV columns and with periodic boundary conditions. At 
the bonds connecting adjacent pair of lattice sites we attach a variable a = 0,1,2,..., 
representing the number of lines (or particles) in the bond. The links (a, (5, 7, 8 = 0, 1,2,.. .) 
connected with a given site produce a vertex configuration, with Boltzmann weight as 
in figure 1. The model we consider obey a type of "ice condition" [1], where the Boltzmann 
weights are chosen such that the number of lines is conserved as we go from the direction 
(a, (3) to (7, 8), i. e., a + [3 = 7 + 8. The non-zero Boltzmann weights are given by 

c.0,0 _ p 0,l _ c.1,0 _ h 
n 0,0 — tt ; -"-0,1 — ^ljO — °5 

R% = Ro°i — c -> (1) 

and 

p0,2 _ pl,l _ pm,n+l _ r p n,0 _ in / n-1 
•"•1,1 — rt 0,2 — •"Ti.m+l — °i -"Vo — u l a 1 

T>m+l,n _ ° 2 (b\ r?m+l,0 _ pi,m _ ( & \ (0 \ 

where m, n > 0, I > and a, b, c are arbitrary parameters. For convenience we use the 
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Baxter's parametrization [3] 

a = sinh -(A — u), b = sinh -(A + u), c = sinhA, (3) 

where u and A are the spectral and crossing parameters, respectively. If we exclude in (l)-(2) 
the configurations where a, (3, 7, 5 > 1 we obtain the standard six- vertex model [3]. 

To calculate the partition function of the system we use the row-to-row transfer matrix 
T, with elements T{ a } { a /}, given by the product of the Boltzmann weights connecting two 
sucessive rows {«i, a 2 , ■ ■ ■ , ocn} and {a[, a' 2 , . . . , a' N } 

00 

T(u) {aUa , } = £ RZ$ N (u)RZ>) ■ --RZt-M ( 4 ) 

/3 1 ,/3 2 ,...,/3 JV =l 

The above transfer matrices with different values of the spectral parameter commute among 
themselves, i. e., 

T{v)T{u) =T{u)T{v). (5) 

A sufficient condition for this commutativity is the existence of (u,u') satisfying the 
Yang-Baxter equations [3] 

E </'W<f {u')X^,{u,u') = £ x(/\u,u')WS / {u')R P 1 iM- (6) 

7,5V 7,<5'V" 

Although we did not succeed in finding X"f(w, w') we have checked equation (5) directly. 
We want to solve the eigenvalue equation of the transfer matrix 

= A^. (7) 

An obvious way to describe the state of a row is to specify the positions of individual particles 
(lines). Let 

$(xi,x 2 ,...,x n ) (8) 

be the amplitude corresponding to the state of a row with particles (lines) at the bonds 
Xi, x 2 , ■ ■ ■ , x n . Unlike the ordinary six- vertex model, some of Xi may coincide 
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xi < x 2 < ■ ■ ■ < x n (9) 

A possible way to represent the model under consideration is to draw as much lines on 
an edge as the value of the variable a (number of particles) of this edge. An example of 
arrangement of lines is shown in figure 2. From the definition of our model (1-2) and Fig. 
2 it is clear that as in the ordinary six-vertex model [0 the lines link together by forming 
continuous noncrossing paths through the lattice. If one starts by following a path upwards 
or to the right, then one will always be travelling in one or other of these two directions, 
never downwards or to the left. The cyclic boundary conditions ensure that a path never 
ends. A typical configuration of lines in adjacent rows is shown in Fig. 3. Consequently 
the total number of particles (lines) on arrow is a conserved quantity and our model can be 
treated similarly as in [1,3]. However in our case the calculation is more complicated. We 
write the amplitude in (8) in the form of the Bethe ansatz [3,4] 

...,x n ) = Y^A-Pt-Pn exp^fcp^i + . • . + ik P „x n ), (10) 
p 

where the sum is over all the permutations P of (1,2, ... ,n) and ki(i = 1,2, ...,n) are 
unknown quasimomenta. The eigenvalues A^} corresponding to the eigenvector ^{k} with 
amplitudes (10) are given by 

n n 

A {k} = a N Y[L(k J ) + b N Y[N(k J ), (11) 
j=i j=i 

where L(k) and N(k) can be calculated by considering a single particle n — 1. In our case 
we have 

L(k) = [ab + (c 2 - b 2 )e ik }/[a 2 - abe l % 

N(k) = e lk [l -c 2 (a 2 -abe ik )- 1 ]. (12) 

Considering the eigenvalue equations at the boundaries of the inequalites (9) after a long, 
but straightforward calculation, we obtain the conditions for the coefficients Af\...p n , 

A... PlP2 ... = -(1 - 2Ae ikp i + e ikp i +ikp z) 
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• (1 - 2Ae ifcp 2 + e ik Pi+ ik p 2 )^e ikp 2- ikp iA...p 2Pl ..., (13) 
A Pl ... Pn = A P2 ,„p nPl exp(ik Pl N), (14) 

where 

A = (a 2 + b 2 - c 2 ) /2ab = - cosh A. 

The Bethe-ansatz equations that come from (13)-(15) and fix the values of {A;} are 
simplified by using the standard parametrization of the six-vertex model [3] 

e i k] = sinh|(A-^) 



sinh | (A + Uj) ' 



(15) 



and are given by 

sinh |(A + Uj] 1 



sinh |(A — u 



3)] 



= (-iy-i e -iP TT smh[\ + l( Uj - ui) 
jLj sinh [A — |(-Uj — u t ) 



(16) 



where 

n 

P = h (mod27r) 
i=i 

is the associated momentum of the eigenstate ^{k}- The eigenvalue (11) is now given by 

n _ h n ~ 

A(«) W = a N II L M + ^(-) n ea ^ P ) II (17) 

where 

sinh[A + \{u — Uj)} 



L( Uj ) 



sinh \ (uj — u) 



~, . sinh[A — ^(u — Uj)] , , 

M(uj) = 18 

sinh ±(u - Uj) 

Equations (10)-(18) completely describe the eigenvectors and eigenvalues of the row-to-row 
transfer matrix of the model with Boltzmann weights (l)-(2). From the comparison of 
these equations with the corresponding equations of the six-vertex model it is clear that the 
eigenvalues of the transfer matrix of the model under consideration up to the multiplicative 



factor a n coincide with the eigenvalues of the transfer matrix of a six-vertex model with 
N + n columns but with Boltzmann weights modified along a vertical " seam" in the lattice 
[13,14]. The boundary angle, that modifies the Boltzmann weights along the seam is equal 
to the momentum P of the state. 

Let us consider the thermodynamic limit where n — > oo and N — > oo, but the density 
of particles (or density of lines) p = n/N is kept fixed. The equations (16)-(17) can be 
rewritten in terms of integral equations. In the region — 1 < A < 1, \ = ip, < p < n we 
have 

/ x f u ° sm(2p)p(u') , . sinu . , , . 

/•M 

/ p{u)du = p, (20) 

J — un 



1 ruo 
— ln[A(u)] = — plna + (1 + p) maxjlna + / \n[L{u)\p{u)du] 
-/V J —uo 

In b + — + / \n[M{u)]p{u)du}. (21) 



Since the maximum eigenvalue has momentum zero it is interesting to compare equations 
(19)-(20) with the corresponding ones in the six- vertex model [1]. This comparison enables 
us to express the free energy f(p) of our model, at a given density p, in terms of the free 
energy f 6v (p) of the six- vertex model, i.e. 

f(p) = (3- 1 p\na+ (1 + p)f^-P-). (22) 

1 + p 

We can also obtain an exactly integrable quantum Hamiltonian, related with our vertex 
model by taking the logarithmic derivative of the row-to-row transfer matrix (4) at u — —A. 
This gives us a new integrable quantum model with Hamiltonian 

L oo 

H = -N cosh A + 2 cosh A £ £ (m - l)Ef' m 

j=l m=l 

L oo 

+ E E (E?' m+1 E%S' n + Ef^E]^ 1 ) (23) 

j=l m,n=l 
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where E™' m are infinite-dimensional matrices with elements given by (E l ' m )ij = 5i^5 m j. This 
Hamiltonian describes the dynamics of particles on a lattice with no volume exclusion. It 
contains a hopping term (second sum) which allows the motion of one of the particles at 
a given site to its neighbouring sites, independently of its relative occupation and a static 
interaction (first sum) proportional to the occupation number at the site, rii = J2m=i mE™' m . 
The line conservation in the vertex model translates into the commutation of (23) with the 
total number of particles n = J2i=i n i- The eigenenergies in a given eigensector with n 
particles are given by E^} = — 2 X^=i cos kj, up to an additive constant which we have 
dropped, where {kj} are the solutions of (16). In the thermodynamic limit the ground-state 
energy is obtained in terms of the solution p(u) of (19)-(20), namely 

Eq n f u ° r sin 2 a , . . , , . 

■77 = 2/ [cosp ^ ]p(u)du. 24 

N J~u l ^ cosh M — cos /i ' v ' 

In figure 3 we show the ground-state energy per site for < p < 1, obtained by solving 
numerically (19), (20) and (24). 

To study the physical properties of this quantum chain we calculated the long-distance 
behaviour of the correlation functions by exploring the consequences of conformal invariance 
(see [15-18] and references therein). The critical fluctuations are governed by a conformal 
field theory with central charge c = 1 and as in the XXZ chain the critical exponents are 
like those of a Gaussian model [19,15]. The long-distance behavior of the density- density 
and of the pair-pair correlation functions are given by 

(p(r)p(O)) ~ p 2 + A ir - 2 + A 2 r- a cos(2irpr) (25) 

and 

(r] + (r)ri(0)) - Br~ p , (26) 

where p(r) = J2m=i m E" l ' m and r](r) = J2m=2 E™ -2 ' 171 are the density and pair operators. 
The exponents a and f3 describing the algebraic decay of the correlations are given by 

a = /3- 1 = 2[eWf, (27) 
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where the dressed charge function is obtained by solving the integral equation 

«°»^L, mS KJ,,T-ci(2^ ^ (1 ^ 1 ' (28) 

As usual in one dimension the model does not show finite off-diagonal long-range order. 
However, the power decay of the pair-pair correlation (27) indicates the preference of pairing 
of particles whenever the exponent (5 of this correlation is smaller than that of the density- 
density correlation a [16]. The exponent (3 is plotted in figure 4 for < p < 1 and some values 
of the parameter //. These results are obtained by solving numerically the integral equations 
(19-20) and (28). We see from this figure that for all values — 1 < A = — cosp < 1 there 
exists a regime with dominant pair correlations. The analogous behaviour of correlation 
functions was also observed in the models considered in Refs. [20-21]. But differently from 
these models, in the present case we have the nice feature that for increasing values of p, the 
exponent (3 for < A < 1 tends toward zero and there is no regime where the density- density 
correlations dominate. 

It should be emphasized that the formalism exposed here for creation of new integrable 
models applies to lattice systems with higher spin like the Perk-Shultz [8] and Sutherland 
models [10] as well as Fateev-Zamolodchikov [11] and Izergin-Korepin models [12]. 

This work was supported in part by CNPq anf FINEP-Brazil. 
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FIGURES 



FIG. 1. The vertex formed by the link variables a, (3, 7, 5 and having Boltmann weight R^'g. 
FIG. 2. An example of a line representation of a configuration in the unlimited vertex model. 

FIG. 3. A typical arrangement of lines of two adjacent rows. 

FIG. 4. The ground-state energy per site Eq/N as a function of the density p = n/N and 
anisotropy — 1 < A = — cos p < 1, with (a) p = n/10, (b) p = 7r/4, (c) p = 7r/3, (d) p = ir/2, (e) 
p = 5tt/6 and (f) p = lOrr/12. 

FIG. 5. The critical exponent (3 as a function of the density p = n/N , and for some values 
of the anisotropy A = —cosp, with (a) p = 7r/10, (b) p = tt/6, (c) p = 7r/4, (d) p = ir/2, (e) 
p = 5vr/6 and (f) p = lOvr/12. 
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